How accurately can the deuterium abundance be deter- 
mined ? 
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I Modern observations of the H+D absorption at high redshift play an 

■ important role in determining the primordial hydrogen isotopic ratio D/H 
= Ndi/Nui (the ratio of the DI to HI column densities). Measurements of 

^ ■ TVdi and TVhi from QSO spectra can provide a sensitive test of the predictions 

I of big bang nucleosynthesis [BBN] if their precision is comparable with the 

■ theoretical uncertainties (~ 15 %). 

I It has been shown, however, that the determination of N depends on the 

. assumptions made with respect to the line broadening mechanism and may 

I be ambiguous especially for the case of optically thick lines [1-8]. The com- 

Q^ ■ monly used procedure is to fit Voigt profiles, accounting for radiative damping 

and Doppler broadening in the microturhulent limit (completely uncorrelated 
^ . bulk motions). But the interpretation of observed profiles may be substan- 

tially changed if a mesoturhulent model is used, i.e. if the influence of a finite 
correlation length in the stochastic velocity field on the line forming process 
^ ' is accounted for. This appears to be appropriate in particular with respect to 

• the H+D absorption observed in QSO spectra. The hydrogen absorbers with 

A^Hi ~ 10^^ — 10^* cm^^ (which are especially useful for deuterium observa- 
tions) are usually associated with extended gas halos of foreground galaxies. 
Direct observations |^ of distant galaxies {z > 2) reveal complex morpholo- 
^ . gies and kinematics of the extended HI Lya absorption gas with projected 

sizes up to 50 kpc. The inferred high values of the rms velocity dispersion 
{at ~ 20 — 50 km s^^) and a few kpc sizes of the absorbing regions make the 
mesoturhulent approach more adequate in this case. 

Within the framework of the mesoturhulent model, one has to determine 
both the physical parameters of the gas and the velocity field structure v{s) 
along the line of sight. This problem may be successfully solved by using a 
Reverse Monte Carlo [RMC] technique 

We applied the RMC procedure to a template H+D Lya profile which 
reproduces the Q 1009+2956 spectrum with the DI Lya line seen at Za = 2.504 
[ p^ . Here we present only a part of our results, the full analysis is given in 

To simulate real data, we added the experimental uncertainties to the 
template intensities which were sampled in equidistant bins as shown in Fig. 1 
by dots and corresponding error bars. A one-component mesoturhulent model 



with a homogeneous density and temperature was adopted. Adequate profile 
fits for three different sets of parameters are shown in panels (al, bl, cl) by 
solid curves, whereas the individual HI and DI profiles are the dashed and 
dotted curves, respectively. The estimated parameters (D/H, TVhi, Tkm, the 
ratio of the rms turbulent velocity to the thermal velocity at/vth, the ratio of 
the cloud thickness to the correlation length L/l) and Xmin values per degree 
of freedom are also listed in these panels for each RMC solution. 

The essential difference between the results of the iwo-component micro- 
turbulent model adopted in |]l0| and ours lies in the estimation of the hydro- 
dynamical velocities in the Za — 2.504 absorbing region. The RMC procedure 
yields at — 25 km s""'^, whereas the iwo-component model leads to crt ~ 2 km 
s~^ which is evidently too low as compared with observations 

In this particular absorption system, the study of the H-l-D Lya profile 
yields D/H = (3.75 ± 0.85) x 10"^ (2(7) which is slightly higher than the value 



2.51^o gg X 10 ^ found in |10 . This difference is, however, very significant 



because it leads to limits on D/H consistent with BBN predictions and obser- 
vational constraints on both extra-galactic ''He mass fraction Yp ||ll| and ^Li 
abundance in the atmospheres of population II (halo) stars ||l^. BBN restricts 
the common interval of the baryon to photon ratio r] for which there is con- 
cordance between the above-mentioned abundances to (3.9 — 5.3) x 10^^° as 
shown by the shaded region in Fig. 2. It implies that 0.014 < ilf,ft,^QQ < 0.020. 

We conclude that the discordance of D /H with the ^He and ^Li measure- 
ments noted in ||l^ is a consequence of the use of the microturbulent model. 
The generalized mesoturbulent model gives good agreement between the light 
element abundances and the BBN predictions. 
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Figure 1: (al, bl, cl) - A template H+D Lya profile (dots with er- 
ror bars) representing the normalized intensities and their uncertainties 
in accord with p^]. The solid curves show the results of the RMC 
minimization, whereas the dotted and dashed curves are the separate 
profiles of DI and HI, respectively. Also shown are the best fitting pa- 
rameters and the values obtained for each case. (a2, b2, c2) - The 
corresponding individual realizations of the velocity distribution u{x) 
in units of at- (a3, b3, c3) - The histograms are the projected velocity 
distributions p{u). 





Figure 2: Comparison of predicted primordial abundances (denoted by 
the subscript p) with observational bounds. The theoretical BBN D, 
^He, and ^Li yields (solid curves) and their uncertainties (dashed curves) 
as function of r] (the baryon to photon number ratio) are computed by 
using the parameterization given in [13|. The abundances of D and ^Li 
are number ratios, whereas Yp is the mass fraction of ^He. The hight of 
the solid-line rectangles in panels b and c give the bounds from recent 
measurements of extra-galactic Yp |11] and ^Li in population II (halo) 
stars [^]. The horizontal widths give the range of ry compatible with 
the measurements. In panel a the solid-line rectangle corresponds to the 
uncertainty region for D/H towards Q1009+2956, as calculated by the 
RMC procedure for an arbitrary velocity field configuration |^]. The 
upper and lower dotted-line rectangles are defined by the confidence 
regions (computed in Q) for the fixed velocity field structures shown 
in Fig. I(c2) and Fig. I(a2), respectively. The shaded region is the 
window for r] common to all results. 



